De novo emergence of a remdesivir resistance mutation during treatment of persistent SARS-CoV-2 infection in an immunocompromised patient: A case report (preprint)

SARS-CoV-2 remdesivir resistance mutations have been generated in vitro but have not been reported in patients receiving treatment with the antiviral agent. We present a case of an immunocompromised patient with acquired B-cell deficiency who developed an indolent, protracted course of SARS-CoV-2 infection. Remdesivir therapy alleviated symptoms and produced a transient virologic response, but her course was complicated by recrudescence of high-grade viral shedding. Whole genome sequencing identified a mutation, E802D, in the nsp12 RNA-dependent RNA polymerase which was not present in pre-treatment specimens. In vitro experiments demonstrated that the mutation conferred a ~6-fold increase in remdesivir IC50 but resulted in a fitness cost in the absence of remdesivir. Sustained clinical and virologic response was achieved after treatment with casirivimab-imdevimab. Although the fitness cost observed in vitro may limit the risk posed by E802D, this case illustrates the importance of monitoring for remdesivir resistance and the potential benefit of combinatorial therapies in immunocompromised patients with SARS-CoV-2 infection.

Clinical performance characteristics of the Swift Normalase Amplicon Panel for sensitive recovery of SARS-CoV-2 genomes (preprint)

Amplicon-based sequencing methods have been central in characterizing the diversity, transmission and evolution of SARS-CoV-2, but need to be rigorously assessed for clinical utility. Here, we validated the Swift Biosciences SARS-CoV-2 Swift Normalase Amplicon Panels using remnant clinical specimens. High quality genomes meeting our established library and sequence quality criteria were recovered from positive specimens with a 95% limit of detection of 40.08 SARS-CoV-2 copies/PCR reaction. Breadth of genome recovery was evaluated across a range of Ct values (11.3 - 36.7, median 21.6). Out of 428 positive samples, 406 (94.9%) generated genomes with < 10% Ns, with a mean genome coverage of 13,545X/SD 8,382X. No genomes were recovered from PCR-negative specimens (n = 30), or from specimens positive for non-SARS-CoV-2 respiratory viruses (n = 20). Compared to whole-genome shotgun metagenomic sequencing (n = 14) or Sanger sequencing for the spike gene (n = 11), pairwise identity between consensus sequences was 100% in all cases, with highly concordant allele frequencies (R2 = 0.99) between Swift and shotgun libraries. When samples from different clades were mixed at varying ratios, expected variants were detected even in 1:99 mixtures. When deployed as a clinical test, 268 tests were performed in the first 23 weeks with a median turnaround time of 11 days, ordered primarily for outbreak investigations and infection control.

A bifluorescent-based assay for the identification of neutralizing antibodies against SARS-CoV-2 variants of concern in vitro and in vivo (preprint)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged at the end of 2019 and has been responsible for the still ongoing coronavirus disease 2019 (COVID-19) pandemic. Prophylactic vaccines have been authorized by the United States (US) Food and Drug Administration (FDA) for the prevention of COVID-19. Identification of SARS-CoV-2 neutralizing antibodies (NAbs) is important to assess vaccine protection efficacy, including their ability to protect against emerging SARS- CoV-2 variants of concern (VoC). Here we report the generation and use of a recombinant (r)SARS-CoV-2 USA/WA1/2020 (WA-1) strain expressing Venus and a rSARS-CoV-2 expressing mCherry and containing mutations K417N, E484K, and N501Y found in the receptor binding domain (RBD) of the spike (S) glycoprotein of the South African (SA) B.1.351 (beta, {beta}) VoC, in bifluorescent-based assays to rapidly and accurately identify human monoclonal antibodies (hMAbs) able to neutralize both viral infections in vitro and in vivo. Importantly, our bifluorescent-based system accurately recapitulated findings observed using individual viruses. Moreover, fluorescent- expressing rSARS-CoV-2 and the parental wild-type (WT) rSARS-CoV-2 WA-1 had similar viral fitness in vitro, as well as similar virulence and pathogenicity in vivo in the K18 human angiotensin converting enzyme 2 (hACE2) transgenic mouse model of SARS-CoV-2 infection. We demonstrate that these new fluorescent-expressing rSARS- CoV-2 can be used in vitro and in vivo to easily identify hMAbs that simultaneously neutralize different SARS-CoV-2 strains, including VoC, for the rapid assessment of vaccine efficacy or the identification of prophylactic and/or therapeutic broadly NAbs for the treatment of SARS-CoV-2 infection.

A SARS-CoV-2 Nucleocapsid Variant that Affects Antigen Test Performance (preprint)

More than one year into a global pandemic, SARS-CoV-2 is now defined by a variety of rapidly evolving variant lineages. Several FDA authorized molecular diagnostic tests have been impacted by viral variation, while no reports of viral variation affecting antigen test performance have occurred to date. While determining the analytical sensitivity of the Quidel Sofia SARS Antigen FIA test (Sofia 2), we uncovered a high viral load specimen that repeatedly tested negative by this antigen test. Whole genome sequencing of the specimen uncovered two mutations, T205I and D399N, present in the nucleocapsid protein of the isolate. All six SARS-CoV-2 positive clinical specimens available in our laboratory with a D399N nucleocapsid mutation and CT < 31 were not detected by the Sofia 2 but detected by the Abbott BinaxNOW COVID-19 Ag Card, while clinical specimens with the T205I mutation were detected by both assays. Testing of recombinant SARS-CoV-2 nucleocapsid with these variants demonstrated an approximate 1000-fold loss in sensitivity for the Quidel Sofia SARS Antigen FIA test associated with the D399N mutation, while the BinaxNOW and Quidel Quickvue SARS Antigen tests were unaffected by the mutation. The D399N nucleocapsid mutation has been relatively uncommon to date, appearing in only 0.02% of genomes worldwide at time of writing. Our results demonstrate how routine pathogen genomics can be integrated into the clinical microbiology laboratory to investigate diagnostic edge cases, as well as the importance of profiling antigenic diversity outside of the spike protein for SARS-CoV-2 diagnostics.

Host-pathogen dynamics in longitudinal clinical specimens from patients with COVID-19 (preprint)

Rapid dissemination of SARS-CoV-2 sequencing data to public repositories has enabled widespread study of viral genomes, but studies of longitudinal specimens from infected persons are relatively limited. Analysis of longitudinal specimens enables understanding of how host immune pressures drive viral evolution in vivo. Here we performed sequencing of 49 longitudinal SARS-CoV-2-positive samples from 20 patients in Washington State collected between March and September of 2020. Viral loads declined over time with an average increase in RT-PCR cycle threshold (Ct) of 0.87 per day. We found that there was negligible change in SARS-CoV-2 consensus sequences over time, but identified a number of nonsynonymous variants at low frequencies across the genome. We observed enrichment for a relatively small number of these variants, all of which are now seen in consensus genomes across the globe at low prevalence. In one patient, we saw rapid emergence of various low-level deletion variants at the N-terminal domain of the spike glycoprotein, some of which have previously been shown to be associated with reduced neutralization potency from sera. In a subset of samples that were sequenced using metagenomic methods, differential gene expression analysis showed a downregulation of cytoskeletal genes that was consistent with a loss of ciliated epithelium during infection and recovery. We also identified co-occurrence of bacterial species in samples from multiple hospitalized individuals. These results demonstrate that the intrahost genetic composition of SARS-CoV-2 is dynamic during the course of COVID-19, and highlight the need for continued surveillance and deep sequencing of minor variants.

Specific allelic discrimination of N501Y and other SARS-CoV-2 mutations by ddPCR detects B.1.1.7 lineage in Washington State (preprint)

Real-time epidemiological tracking of variants of interest can help limit the spread of more contagious forms of SARS-CoV-2, such as those containing the N501Y mutation. Typically, genetic sequencing is required to be able to track variants of interest in real-time. However, sequencing can take time and may not be accessible in all laboratories. Genotyping by RT-ddPCR offers an alternative to sequencing to rapidly detect variants of concern through discrimination of specific mutations such as N501Y that is associated with increased transmissibility. Here we describe the first cases of the B.1.1.7 lineage of SARS-CoV-2 detected in Washington State by using a combination of RT-PCR, RT-ddPCR, and next-generation sequencing. We screened 1,035 samples positive for SARS-CoV-2 by our CDC-based laboratory developed assay using ThermoFishers multiplex RT-PCR COVID-19 assay over four weeks from late December 2020 to early January 2021. S gene dropout candidates were subsequently assayed by RT-ddPCR to confirm four mutations within the S gene associated with the B.1.1.7 lineage: a deletion at amino acid (AA) 69-70 (ACATGT), deletion at AA 145, (TTA), N501Y mutation (TAT), and S982A mutation (GCA). All four targets were detected in two specimens, and follow-up sequencing revealed a total of 10 mutations in the S gene and phylogenetic clustering within the B.1.1.7 lineage. As variants of concern become increasingly prevalent, molecular diagnostic tools like RT-ddPCR can be utilized to quickly, accurately, and sensitively distinguish more contagious lineages of SARS-CoV-2.